Two-dimensional anisotropic Heisenberg antiferromagnet in a field 
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The classical, square lattice, uniaxially anisotropic Heisenberg antiferromagnet in a magnetic 
field parallel to the easy axis is studied using Monte Carlo techniques. The model displays a long- 
range ordered antiferromagnetic, an algebraically ordered spin-flop, and a paramagnetic phase. The 
simulations indicate that a narrow disordered phase intervenes between the ordered phases down to 
quite low temperatures. Results are compared to previous, partially conflicting findings on related 
classical models as well as the quantum variant with spin S=l/2. 

PACS numbers: 75.10.Hk, 75.40.Mg, 75.40.Cx, 05.10.Ln 



I. INTRODUCTION 

The square lattice, uniaxially anisotropic Heisenberg 
antiferromagnet in an external field H parallel to the 
easy axis has been studied since more than two decades 
both in its classical versioni^S*^ and in the quantum 
variant. 5,6 However, important features of the phase di- 
agram have not been definitely clarified. 

The anisotropy may be introduced in various ways. In 
the XXZ case, the Hamiltonian may be written as 
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SfSfi-H^S*, (1) 



where for the classical model Sf,Sf,S? are the compo- 
nents of a unit vector corresponding to the spin at site 
i of the lattice, and the sum runs over all nearest- 
neighbor pairs. The coupling constant J and the field 
H are positive; the anisotropy parameter A may range 
from zero to one. The model is known to exhibit, for 
< A < 1, both for simple cubic and square lattices, an 
antiferromagnetic (AF) phase at low temperatures and 
low fields, a spin-flop (SF) phase with spins canted to- 
wards the z-axis at larger fields, and the disordered, para- 
magnetic phase. 

For the simple cubic lattice, the three phases are be- 
lieved to meet at a bicritical point belonging to the 
universality class of the isotropic Heisenberg model with 
0(3) symmetry. Elsewhere, the boundaries to the disor- 
dered phase display either Ising criticality in the AF case 
or XY criticality in the SF case. 

For the square lattice, different scenarios on the fate 
of the bicritical point have been discussed following the 
early analysis of Landau and Binderi In particular, (a) 
the bicritical point may move to zero temperature, in ac- 
cordance with the well-known theorem by Mermin and 
Wagner, 8 with a presumably narrow disordered phase in 
between the AF and SF phases; (b) at low temperatures, 
there may be a direct transition of first order between 
the AF and SF phases, with a tricritical point on the 
boundary of the AF to the paramagnetic phase; and (c) 
a 'biconical phase', in which the spins on only one of 
the sublattices of the antiferromagnet are canted, may 
intervene between the AF and SF phases at low tem- 
peratures. However, based on the early simulations^ at 



A = 4/5, none of these three scenarios has been strongly 
or even definitely favored. 

More recently, related two-dimensional antiferromag- 
nets in a field have been studied. For a classical antiferro- 
magnet with nearest-neighbor interactions and a single- 
ion anisotropy it has been suggested that the bicritical 
point occurs at zero temperature. However, evidence is 
provided by analytic approximations which are in quanti- 
tatively rather poor agreement with Monte Carlo data. 4 
A phase diagram with a tricritical point and a direct 
transition between the AF and SF phases has been deter- 
mined for the quantum variant of the XXZ-model with 
spin S = 1/2 and A = 2/3, see Eq. Q, correspond- 
ing to the hard-core boson Hubbard models For an ex- 
perimentally motivated classical model 9 with single-ion 
anisotropy and further neighbor couplings a topologically 
similar phase diagram has been obtained^ 

Experimentally, there seem to be quite a few quasi two- 
dimensional antiferromagnets with uniaxial anisotropy, 
including, for instance, Rb2MnF4, Rb2MnCi4, K 2 MnF4, 
La 5 Ca 9 Cu2404i, and Mn(HCOO)2. 2 i 3 ' 9 i u ' 12 i 13 ' 14 How- 
ever, the above sketched subtleties of the phase dia- 
gram have not been fully elucidated, perhaps due to 
additional interactions like interlayer couplings, or addi- 
tional anisotropies like the breaking of the XY isotropy, 
which, even when being weak, are expected to affect 
significantly critical properties^ Inevitable defects may 
eventually play an important role in two-dimensional an- 
tiferromagnets in a field being related to random-field 
systems li 5 - 

In this article, we shall present results of extensive 
Monte Carlo simulations on the classical XXZ Heisen- 
berg antiferromagnet on a square lattice, mostly setting 
the anisotropy parameter to A = 4/5 and A = 2/3. Our 
findings indicate that a narrow disordered phase inter- 
venes between the AF and SF phases down to quite low 
temperatures. 

The layout of the article is as follows. In the next 
section, basic properties of the model will be discussed 
and the quantities computed in the simulations will be 
listed. The transition from the SF to the paramagnetic 
phase will be discussed in Sec. IIIII followed by a section 
dealing with the boundary line of the AF phase. A brief 
discussion and summary will conclude the paper. 
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II. BASIC PROPERTIES AND QUANTITIES OF 
INTEREST 

At zero temperature, T = 0, the XXZ model on 
a square lattice, see Hamiltonian (1), may be easily 
solved exactlyi The antiferromagnetic structure, being 
the ground state at low fields, becomes unstable against 
the spin-flop state at 

H cl = 4Ja/1- A 2 . (2) 

Increasing the field, the paramagnetic state gets stable 
at H > H C 2, with 

ff c2 = 4J(l + A). (3) 

The tilt angle, 9, between the z-axis and the direction 
of the spins in the SF state, H c \ < H < H C 2, is given by 

9 = cos- 1 [ J ff/4J(l + A)]. (4) 

Obviously, in the isotropic limit, A = 1, one gets 
H c i = and H C 2 = 8 J, with the Ising-like AF state be- 
ing squeezed out, the ground state at H = having the 
rotationally invariant 0(3) symmetry. In the Ising limit, 
A = 0, one has ti c \ = H C 2 = 4J, i.e. the SF structure is 
certainly not a ground state. 

At low temperatures, considering < A < 1, the AF 
and SF states give rise to ordered phases. In the AF 
case, one expects long-range order with the z-component 
of the sublattice magnetization as the order parameter. 
The transition to the paramagnetic phase is believed to 
be continuous and of Ising type, at least at small fields. 
In the SF phase, the rotational invariance of the x- and 
y-components of the canted spins may lead to algebraic 
order and a transition of Kosterlitz-Thouless type to the 
disordered phaser^i^ as it is known to hold, for instance, 
for the two-dimensional XY model. 

Phase diagrams of the XXZ model, as obtained from 
our Monte Carlo simulations, are depicted in Figs. 
andEJ}. We set A = 4/5 and A = 2/3, to allow for 
comparison with previous findings^fi The diagrams will 
be discussed in detail later. Unless indicated otherwise, 
here and in the following error bars in the figures are 
smaller than the size of the symbols. 

In the Monte Carlo (MC) simulations, we consider 
square lattices with L x L sites, employing full periodic 
boundary conditions. The linear dimension L ranges 
from 2 to 240, especially to study finite-size effects al- 
lowing for extrapolation to the thermodynamic limit, 
L — > oo (see below). Applying the standard Metropolis 
algorithm, each MC run consists of at least 10 6 (and up 
to 10 8 ) Monte Carlo steps per site. To obtain averages 
and error bars, we usually take into account several, typ- 
ically about ten (and up to forty), realizations choosing 
various random numbers. In selected cases, especially at 
low temperatures close to the boundaries of the AF and 
SF phases, the MC runs are started with different initial 
configurations to check for correct equilibration. 
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FIG. 1: Simulated phase diagrams of the XXZ model on a 
square lattice, for the values (a) A = 2/3, and (b) A = 4/5 
of the anisotropy parameter. Squares denote transitions be- 
tween the SF and paramagnetic phases, circles refer to the 
boundary of the AF phase. Lines, here and in the following 
figures, are guides to the eye. 

We compute quantities of direct experimental inter- 
est as well as other quantities which enable us to con- 
veniently determine the phase transition lines and crit- 
ical properties. In particular, we calculated the specific 
heat C, both from the energy fluctuations and from the 
temperature derivative of the energy. Various magneti- 
zations were computed: Especially, we recorded (i) the 
^-component of the total magnetization, 

<M Z >=(J2S1)/L 2 ; (5) 

i 

(ii) the square of the ^-component of the staggered mag- 
netization < (M|) 2 > (or, similarly, the absolute value) 
to describe the order in the antiferromagnetic phase, 

< (M*) 2 >= ([XX -E^] 2 )A 2 ' (6) 

i a ib 
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summing over all sites, i a and ib, of the two sublat- 
tices a and b of the antiferromagnet; as well as (iii) the 
square of the staggered transverse sublattice magneti- 
zation < (Mf u ) 2 >, to describe the ordering in the SF 
phase, 

< (Mf y ) 2 > = < (A/f) 2 > + < (Mf) 2 >, (7) 

in full analogy to Eq. © . Alternately, one may compute 
the sum of the squares of each sublattice magnetization 
for both transverse components, as before^ 

In addition, we monitored the magnetic (staggered) 
susceptibility X( s ) i which may be obtained from the fluc- 
tuations of the (staggered) magnetization or its deriva- 
tive with respect to the (staggered) field. To identify the 
type of transition along the boundary of the AF phase, 
the fourth-order, size dependent cumulant U[ of the stag- 
gered magnetization, the Binder cumulant^i is supposed 
to be rather useful: 



1- < (M*) 4 > L I (3 < {Mif >1) 



(8) 



where < (Mf) 4 > is defined in analogy to < (M|) 2 >. 
The corresponding cumulant U^ v may be useful to study 
the boundary of the SF phase. To identify the type of the 
phase transition, we also study the fourth-order cumulant 
of the energy^ 



V L = 1- < £ 4 >l / (3 < E 2 >\) 



(9) 



with E being the energy per site. 

Further relevant information on the phase transitions 
may be inferred from histograms of the order parameter 
and of the energy. 



III. TRANSITION FROM THE SPIN-FLOP TO 
THE PARAMAGNETIC PHASE 



For square lattices, the spin-flop phase has been ar- 
gued to be of Kosterlitz-Thouless type; 1 ' 4 ! 6 where trans- 
verse spin correlations, i.e. < SfSf, + SfSf, >, decay 
algebraically with distance \i — i'\ for widely separated 
sites i and i' . Closely related, the transverse sublattice 
magnetization < (M^ v ) 2 >, see Eq. Q, describing the 
ordering in the SF phase, is expected to behave for T > 
and sufficiently large systems as 



(10) 



with g approaching 1/4 at the transition from the SF 
to the paramagnetic phasejiS* 1 ^ and g = 2 in the para- 
magnetic phase. Thence, in two dimensions, the order 
parameter, < (M^ v ) 2 >, vanishes in the SF phase as 
L — > oo at all temperatures T > 0. 

In fact, as illustrated in Fig. [3 the MC data clearly 
show that the magnetization decays with system size 
both in the SF and the paramagnetic phase. To de- 
termine the boundary of the SF phase, Tsf, from the 



size dependence of the transverse sublattice magnetiza- 
tion, Eq. (10), one may study the effective exponent, as 
usual 3 



g e s(L) = - 



din < (M*yy > L 

dlnL 



(11) 



in its discretized form, comparing data for consecutive 
system sizes, L\ and L 2 , L 2 > L\, with 

(1 s ln[< {M-v) 2 > L2 I < (MP) 2 > Ll ] 

.9cff £o = VTTTT\ ' ( 12 ) 

ln(L 2 /£i) 



where Lq = \JL\L2. Indeed, when crossing the phase 
boundary by, for instance, fixing the field and decreasing 
the temperature (see Figs, la and lb), for large systems 
<7eff tends to drop rapidly from rather high values ap- 
proaching 2, characterizing the decay in the disordered 
phase, to a small value close to 1/4 at the transition to 
the SF phase. Deeper in the SF phase, g e g decreases 
somewhat to even lower values. 
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FIG. 2: Square of the staggered transverse sublattice magne- 
tization, < (Mf y ) 2 >l, vs. temperature at fixed field H/J = 
6.0 for systems of size L = 10, 20, 40, 80, and 120 (from 
top to bottom), close to the boundary between spin-flop and 
disordered phases. The anisotropy parameter is A = 4/5. 

This behavior of the effective exponent is exemplified 
in Fig. |21 allowing one to estimate Tsf- We checked the 
approach by analyzing the two-dimensional (planar) XY 
model, reproducing quite accurately the transition tem- 
perature as obtained from elaborate numerical studies^i 

The Kosterlitz-Thouless character of the transition be- 
tween the SF and the paramagnetic phases is also re- 
flected in the thermal behavior of the specific heat C, 
which displays a non-critical maximum close to, but not 
exactly at the transition. Of course, from simulational 
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data one cannot identify the expected essential singular- 
ity of C at the transition. 



IV. THE BOUNDARY LINE OF THE 
ANTIFERROMAGNETIC PHASE 
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FIG. 3: Effective exponent <? e ff, see Eqs. (ff) and (f2), vs. 
temperature at fixed field H/J — 2.7, comparing systems of 
size (Ln,Li)= (20,10), (30,20), (40,30), (60,40), and (80,60) 
(from bottom to top), close to the boundary between spin-flop 
and disordered phases. The anisotropy parameter is A = 4/5. 

Approximate analytic expressions for the boundary be- 
tween the SF and the paramagnetic phases may be ob- 
tained by the following considerations, slightly modifying 
previous arguments. 1 Using polar coordinates and fixing 
the z-component of the spins in the SF phase to its value 
in the ground state, Eq. J3J, the transition temperature 
Tsf may be approximated bj*i 



ksTsF k B T KT 



J 



J 



A [1 - {H/H c2 f} 



(13) 



where 1cbTkt/J refers to the Kosterlitz-Thouless transi- 
tion temperature of the XY model with coupling J, i.e. 
the classical isotropic model with a two-component spin 
vector of length one^i Indeed, this approximate expres- 
sion may be regarded as an upper bound to the true 
transition temperature, because the fluctuations in the 
z-component (which, in turn, are coupled to the fluctua- 
tions in the transverse spin components) tend to reduce 
the transition temperature. In addition, the change of 
the tilt angle with temperature, at fixed field, will 
affect the transition temperature. In fact, Tgp, as de- 
termined in the simulations, see Figs.^i and^ 3 , is seen 
to be lower than suggested by Eq. ifHfl) . In any event, 
the approximation does not hold in the vicinity of H&, 
where the fluctuations of the longitudinal and transverse 
spin components are strongly correlated. Among others, 
that region will be discussed in the next section. 



We determined the boundary line of the AF phase 
by monitoring especially the specific heat C, the square 
of the z-component of the staggered magnetization, 
< (Mg ) 2 >, the staggered susceptibility, xfi an d the 
Binder cumulant, £/£. The transition temperatures fol- 
low from finite-size extrapolations of the Monte Carlo 
data. Results for A = 4/5 and A — 2/3 are displayed in 
Fig.E 

At low fields and high temperatures the transition is 
expected to be continuous and of Ising type. Indeed, 
we observe, for instance, a logarithmic divergence in the 
height of the peak in the specific heat as a function of 
system size, C max oc In L, and an effective critical ex- 
ponent of the order parameter /3 e g consistent with the 
asymptotic value 1/8 when approaching the transition, 
fixing the field and varying temperature. 

The transition line of the AF phase shows a maximum 
in the critical field, as a function of temperature, being 
somewhat higher than H c i, Eq. both for A = 2/3 
and A = 4/5 (a feature which has not been mentioned 
in the early worhi). In that part of the phase diagram, 
see Fig. 0] for A = 4/5, the boundary lines of the SF 
and AF phases approach each other, when lowering the 
temperature. 
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FIG. 4: Details of the phase diagram of the XXZ model on a 
square lattice with A = 4/5. Squares refer to the boundary 
of the SF, circles to that of the AF phase. 

Fixing the field to be slightly below its maximal crit- 
ical value and decreasing the temperature, one observes 
a rather asymmetric behavior in the staggered magne- 
tization < (M*) 2 >, as illustrated in Fig. The order 
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parameter of the AF phase rises fairly gradually enter- 
ing the AF phase from the paramagnetic phase, while it 
drops down rather rapidly on approach to the SF phase. 
This asymmetry signals either a change in the type of 
the transition at the lower temperature, becoming, pos- 
sibly, of first order, or in the extent of the critical region 
becoming, possibly, very narrow. 




FIG. 5: Staggered magnetization < (M*) 2 > vs. temperature 
at field H/J = 3.0 for A = 2/3 with L = 128 (squares), and 
at H/J= 2.42 for A = 4/5 with L = 120 (circles), close to 
the maximum of the boundary of the AF phase. 

To study the transition along the boundary of the AF 
phase in more detail, the Binder cumulant £/£ may be 
quite useful, as had been employed beforeAi^ In the 
thermodynamic limit the value of the cumulant at the 
transition point, the critical cumulant U[ =00 (Taf), is 
believed to characterize the type and universality class 
of the transition. From simulational data, the critical 
cumulant may be estimated from the intersection value 
U 1 (L) of the cumulant for different system sizes L and 
L' = bL (in the following, we set b — 2), doing finite- 
size extrapolations i± In Fig. we depict results for 
A = 4/5, monitoring U 1 (L) as a function of tempera- 
ture in the vicinity of the AF boundary line. The inter- 
section value U 1 seems to display one plateau at higher 
temperatures and another plateau, with a smaller height, 
at lower temperatures, with a fairly rapid change in be- 
tween. The plateaus tend to get more pronounced and 
the change becomes sharper as the system size L is in- 
creased. In fact, the plateau at high temperatures is ob- 
viously related to the critical cumulant in the univer- 
sality class of the two-dimensional Ising model, where 
U l=00 (Taf) ~ 0.6106^ Note that U[, computed ex- 
actly at the transition (Taf, -Haf), approaches this value 
closely with only minor finite-size corrections for low 



fields and high temperatures. However, near the max- 
imum of the boundary line of the AF phase, significant 
corrections are observed, which have not been discussed 
in the related quantum variant^ The critical cumulant 
of the apparent plateau at low temperatures tends to in- 
crease weakly with T, being close to 0.4 at A = 4/5, 
for the system sizes we considered, L < 120. Perhaps 
interestingly, at the critical field H c i(A), the cumulant 
seems to approach, at least for small lattices, in the limit 
T — ► 0, a value close to 0.38. With decreasing A, this 
limiting value decreases. 
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FIG. 6: Binder cumulant (7 / (L), as obtained from the in- 
tersection points of cumulants for lattices of sizes L and 2L 
close to the boundary of the AF phase, as a function of tem- 
perature, with L= 10, 20, and 40 (from bottom to top). The 
anisotropy parameter is A = 4/5. 

Note that the turning point of the intersection value 
U 1 (L) is shifted towards lower temperatures, T tu , as L 
increases, see Figs. and \7\ Obviously, there are un- 
certainties in extrapolating those data to the thermody- 
namic limit. A naive linear extrapolation yields an es- 
timate of a temperature, T* , of about 0.4 (in units of 
J/ks), at which non-Ising criticality may set in, possibly 
due to a transition of first order. 

This interpretation, however, has to be viewed with 
much care. Indeed, the fourth order energy cumulant 18 
gives no hint of a transition of first order at the AF 
boundary, at quite low temperatures, 0.3 J/ Ub < T < T* . 
Its minimum near the transition, becoming even more 
shallow at lower temperatures, tends to vanish for larger 
lattices, as expected for a continuous transition. 

Even more strikingly, the size dependences of the spe- 
cific heat as well as of the staggered susceptibility sug- 
gest that the boundary line of the AF phase is still in 
the Ising universality class well below the anticipated 
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FIG. 7: Size dependence of the turning point, T tu , in the 
Binder cumulant U 1 {Li) comparing lattices of linear dimension 
L and 2L, for A = 4/5, see Fig. E] with L ranging from 4 to 
60. 



change at T* . This behavior is exemplified in Figs. |S| 
and El for the case A = 4/5 and H/J = 2.42, where 
ksTAF/J ~ 0.33, i.e. in the part of the phase diagram 
where the AF and SF phase boundaries are hardly dis- 
cernible, see Fig. 0] We find consistency with Ising criti- 
cality, especially, C max oc InL, and the peak of the stag- 
gered susceptibility, (xf)max, grows like L 1 ^ — L 7 / 4 . We 
checked equilibration of the MC data by monitoring time 
series of the energy and the staggered magnetization as 
well as by calculating the specific heat from the energy 
fluctuations and the temperature derivative of the en- 
ergy. The findings suggest that there is a very narrow 
paramagnetic phase intervening between the AF and SF 
phases at least down to that temperature, fc^T/ J w 0.33. 
At even lower temperatures, this type of analysis would 
require very large computational resources to obtain re- 
liable MC data. 

For A = 2/3, similar features are observed in the cor- 
responding part of the phase diagram, where, e.g., a loga- 
rithmic size dependence of the height of the specific heat 
at temperatures close to /cbTaf/^ ~ 0.28 is found, com- 
pare Fig. ^ 

In contrast, results of our simulations of the XXZ 
model on a simple cubic lattice, setting A = 4/5, show 
at temperatures below the bicritical point a direct tran- 
sition of first order between the AF and SF phases, as 
expectediiSi For instance, the maximum of the specific 
heat increases with a power law in L, with an effective 
exponent, for moderate system sizes, 8 < L < 20, being 
significantly larger than the value characteristic for the 
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FIG. 8: Maximum in the specific heat vs. logarithm of the sys- 
tem size L at H/J — 2.42 and A = 4/5 at temperatures close 
to k bT af / J ~ 0.33, see Fig. 4. The dashed line indicates a 
logarithmic size dependence. 

universality class of the three-dimensional Ising model. 
The behavior is, indeed, indicative of a first order tran- 
sition, where C max oc L 3 . Likewise, energy histograms 
show the usual properties of a transition of first order, 
with an overlap of two distinct Gaussian peaks. The 
critical energy cumulant in the three-dimensional case is 
observed to be close to, but distinct from 2/3. 

Note that at fixed low temperature and varying the 
field, the z-component of the magnetization changes 
rapidly in both sublattices simultaneously, within our 
field resolution, both for A = 4/5 and A = 2/3, con- 
tradicting a biconical intermediate phase with a canting 
of the spins in only one sublattice. 

V. DISCUSSION AND SUMMARY 

We have studied the classical, square lattice, uniaxially 
anisotropic Heisenberg antiferromagnet, the XXZ model, 
in a magnetic field, doing extensive Monte Carlo simu- 
lations. We mainly considered two cases of fairly weak 
anisotropy, A = 4/5 and A = 2/3. 

The model displays an antiferromagnetic, a spin-flop, 
and a paramagnetic phase. 

The transition from the antiferromagnetic to the para- 
magnetic phase is seen to be of Ising type at small fields, 
remaining of that type at higher fields and rather low 
temperatures (down to at least ksT/J w 0.33 in the 
case A = 4/5), as inferred, especially, from the size de- 
pendence of the maxima in the staggered susceptibility 
and the specific heat. A naive analysis of MC data of 
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FIG. 9: Doubly logarithmic plot of the size dependence of the 
maximum in the staggered susceptibility at H/J = 2.42 and 
A = 4/5. The dashed line refers to a slope with y/v = 7/4. 

the Binder cumulant along the boundary of the antifer- 
romagnetic phase may lead to a different, but supposedly 
erroneous conclusion. Presumably, the critical region at 
the AF phase boundary becomes rather tiny when the SF 
and AF phase boundaries approach each other, demand- 
ing very large lattices to identify the universality class 
correctly, at least when studying the Binder cumulant. 

The transition between the spin-flop and the disor- 
dered phase is of Kosterlitz-Thouless type. It may be 
located accurately by analyzing the finite-size behavior 
of the transverse staggered magnetization, vanishing with 
increasing system size with a characteristic power law at 
the transition. 

In the quantum variant of the model with spin S = 1/2 
a tricritical point on the AF boundary and a direct tran- 
sition between the AF and SF phases at low tempera- 
tures have been obtained^ Near the tricritical point the 



AF and SF phases are still well separated. At present, 
we can only speculate whether this, perhaps, surprising 
discrepancy between the classical and quantum models 
could be resolved when the classical version would be 
analyzed at even lower temperatures or when MC data 
on the quantum version would be reanalyzed. 

In the isotropic limit, A = 1, one finds, in non- 
vanishing field, H > 0, a spin-flop phase of Kosterlitz- 
Thouless character both in the classical and in the quan- 
tum modeliSi For A < 1 and H = 0, an Ising type 
transition occurs at a non-zero temperature, with the 
transition temperature moving to zero as the anisotropy 
vanishes, A — ► 1, again both in the classical and the 
quantum ca.se 

One of the remaining open crucial questions is, whether 
there is a typical phase diagram for square lattice, uniax- 
ially anisotropic antiferromagnets. In the case of the an- 
tiferromagnetic nearest-neighbor Hcisenberg model with 
single-ion anisotropy, a phase diagram with a paramag- 
netic phase between the antiferromagnetic and spin-flop 
phases extending down to arbitrarily low temperatures 
has been suggested^ In another study of that model, a 
conflicting scenario with a direct transition between the 
AF and SF phases has been favored, presenting, how- 
ever, only few Monte Carlo data, 3 Including more than 
nearest-neighbor interactions in an antiferromagnet with 
single-ion anisotropyjSiifi or taking into account quantum 
fluctuations^ phase diagrams with a tricritical point and 
a direct transition between the AF and SF phases have 
been obtained. 

In any event, we should like to encourage future work 
on this and related models, which may also serve as a 
guide in interpreting experiments on corresponding quasi 
two-dimensional anisotropic antiferromagnets. 
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